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We have studied D^-k^ and — + DfK^ decays using 23.6 fb~^ of data collected at the 

T{5S) resonance with the Belle detector at the KEKB e"^e~ collider. This highly pure — > D~n'^ 
sample is used to measure the branching fraction, B{B'^ -D^tt^) ~ [3.67lo 33(stat.)lQ42(syst.) ± 
0.49(/s)] X 10"^ [fs = Ngi,)-i,)/N^^) and the fractions of B° event types at the T{5S) en- 
ergy, in particular Ng,g,/N^(,)git) = (9O.II4 ± 0.2) %. We also determine the masses 

M{B°) = (5364.4 ± 1.3 ±0.7) MeV/c^ and M{B*) = (5416.4 ± 0.4 ± 0.5) MeV/c^ In addition, 
we observe B^—>DfK^ decays with a significance of 3.5 cr and measure B{B^^~*DfK^) = 
[2.4ll g(stat.)±0.3(syst.) ± 0.3(/s)] x lO"*. 

PACS numbers; 13.25.Hw, 13.25.Gv, 14.40.Gx, 14.40.Nd 



The decay D^tt^ \l] has a relatively large 

branching fraction and is a primary normalization mode 
at hadron colliders, where the absolute production rate 
of B'^ mesons is difficult to measure directly. It pro- 
ceeds dominantly via a Cabibbo-favoured tree process. 
The decay B'^ D~7r+ proceeds through the same tree 
process but may also have additional contributions from 
M^-exchange, so a comparison of the partial widths of 
the two decays can give insight into the poorly known 
T^-exchange process. The Cabibbo-suppressed mode 
DfK"^ is mediated by 6 ^ c and b ^ u tree tran- 
sitions of similar order , in the Wolfenstein param- 
eterization 2]), which raises the possibility of measuring 
time-dependent CP- violating effects It has recently 
become possible to produce 5° events from e+e" colli- 
sions at the T(5S) resonance in sufficiently large numbers 
to achieve interesting and competitive measurements. 
T(5S) events may also be used to determine precisely the 
masses of B* and _B^; the mass difference can be com- 
pared with that of B*" and i?° to test heavy-quark sym- 
metry Q, which predicts equality between them. Proper- 
ties of the T(55) such as the fraction of events containing 
a _B° and the relative proportions of B'^B^, B*B2, and 
BtB* provide additional tests of heavy quark theories 

m. 

In this Letter, we report measurements performed 
with fully reconstructed i?° — > DJtt^ and — > DfK^ 
decays in Lint = (23.6 ± 0.3) fb^^ of data collected 
with the Belle detector at the KEKB asymmetric-energy 
(3.6 GeV on 8.2 GeV) e+e" coUider operated at 
the T(5S') resonance. The beam energy in the 

center-of-mass (CM) frame is measured to be = 
Vs/2 = 5433.5 ± 0.5 MeV with T(5S') ^ TjlS") 7r+ tt", 
T(IS') (Lt^M^ decays 0. The total bb cross sec- 



tion at the T{5S) energy has been measured to be 
^T(5S) ^ (0.302 ±0.014) nb 0, which includes 5°, 
_B+ and _B° events. Three 5° production modes are 
kinematically allowed: B*B^, and B*B*. The 

B* decays electromagnetically to B^, emitting a photon 
with energy i?^ ~ 53 MeV. The fraction of bb events 
containing a Bi*'' Bs*'' pair has been measured to be 
fs = iV^(.)^(.)/7Vfc5 - (19.5t^-0)% g. The number of B° 

mesons in the sample is thus Nga — 2 x Lint x '''^5'^^'' ^ 
fs = {2.78to.tl) X 10^. The 5° production mode ra- 
tios are defined as fs'S* = ^B'^b*/ ^si^^Bi*^' ^B'BO = 
^B*s?/^s<*)_B<*' ^^'^ ^B°B° = Ngogo/N^Mgi,). Belle 
previously measured /b'S" — (93lg) % [ij. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a cen- 
tral drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of 
time-of- flight scintillation counters (TOF), and an elec- 
tromagnetic calorimeter comprised of CsI(Tl) crystals lo- 
cated inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to detect and to 
identify muons. The detector is described in detail else- 
where (ill] . 

Reconstructed charged tracks are required to have a 
maximum impact parameter with respect to the nom- 
inal interaction point of 0.5 cm in the radial direction 
and 3 cm in the beam-axis direction. A likelihood ratio 
Uk/^ = Ck/{C^ + Ck) is built using ACC, TOF and 
CDC [dE/dx) measurements. A track is identified as a 
pion if TZk/t! < 0.6 or as a kaon otherwise. With this 
selection, the identification efficiency for pions (kaons) is 
about 91% (85%), while the fake rate is about 9% (14%). 
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Neutral kaons are reconstructed via the decay Kg 
7r~'"7r~ with no identification requirements for the two 
charged pions. The Kg candidates are required to have 
an invariant mass within ±7.5 MeV/c^ (i4cr) of the 
nominal Kg mass (all nominal mass values are taken 



from Ref. I4I). Requirements on the Kg vertex dis- 



placement from the interaction point and on the dif- 
ference between vertex and Kg flight directions are ap- 
plied. The criteria are described in detail elsewhere [l3(. 
The K*'^ {(f)) candidates are reconstructed via the decay 
K*^ K~^TT~ {4> K^K~) with an invariant mass 
within ±50 MeV/c^ (±12 MeV/c^) of the nominal mass. 

Candidates for arc reconstructed in the three 
modes ^tt", Dj K*°K~-, and Dj K^K- 

and required to have mass within ±15 MeV/c^ (iSa) 
of the nominal Dj mass for D~tt^ and within 

±8 MeV/c^ for B° DfK^. Following Ref. \m], the 
signals for DjTr+ and B° DfK^ are observed 

using two variables: the beam-constrained mass of the 

B^ candidate Mbc = \/-^b'^ ~ Pb" ^^'i t^'^ energy dif- 
ference Ai? = Ego - i?b 1 where (£'^0 , P^o ) is the four- 
momentum of the i?s candidate expressed in the CM 
frame. We select candidates with Mbc > 5.3 GeV/c^ 
and -0.3 GeV < A£; < 0.4 GeV. In each event the B° 
candidate with the Dj mass closest to its nominal value 
is selected for further analysis; only w 1% of events have 
more than one candidate. 

Further selection criteria are developed using Monte 
Carlo (MC) samples based on EvtGen [li] and GEANT 
[ist detector simulation. The most significant source of 



background is continuum events, e+e 



uu, dd, ss, cc. 



In addition, for the DfK^ mode there is also 

a large background from i?" D^tt^, where the 



TT 



+ 



is misidentified as a K^. The expected con- 
tinuum background, A'^bkg, is estimated using MC- 
generated continuum events representing three times the 
data. The expected signal, A'sig, is obtained assuming 
B {B° DJTI+) = 3.0 X 10-3 and fB*Bt = 93% for the 
B° ^ analysis and B {B° DfK"^) = 3.7x10-* 

for the B° DfK^ analysis. For DJK^, we as- 

sume the values of B{B^g D^tt^) and /b-s* obtained 
in the B^ Djir^ analysis. 

To improve signal relative to background, criteria are 
chosen to maximize Nsi^/ Ns\ 



bkg 



evaluated in 

the B*B* signal region (Fig. [1]). Two topological vari- 
ables are used. First, we use the ratio of the second 
and zeroth Fox- Wolfram moments |ljQ|, R2, which has 
a broad distribution between zero and one for jet-like 
continuum events and is concentrated in the range be- 
low 0.5 for the more spherical signal events. Candidates 
for B^ D^TT^ (_B^ DfK^) are required to have 
i?2 < 0.5 (< 0.4). We then use the hehcity angle ^hci 
of the (pTT^ {D^ K*^K^) decays, defined as 

the angle between the momentum of the positive daugh- 
ter of the (/) [K*^) and the momentum of the in the 



{K*") rest frame; for signal decays consisting in a spin- 
particle decaying into a spin-1 particle and a spin-0 
particle, the distribution is oc cos^ 0hoi, while for combi- 
natorial background under Ds signal it is flat. Candi- 
dates for Dj (f>7T~ and — + K*'^K~ are required 
to satisfy Icoseihoil > 0.2 (> 0.35) for the B° D-Tr+ 
{B^ DfK"^) mode. These two selections reject 43% 
(73%) of the continuum while retaining 95% (85%) of the 
B° D-TT+ {B° DJK^) signal. MC studies show 
that background from B'^ and B^ decays is small and 
flat enough to be described together with the continuum 
events for the i?^ — ^ D^n^ mode and is negligible for 
the — > DJK^ mode. The most relevant background 
from decays is i?° D*~tt+ . 

For each mode, a two-dimensional unbinned extended 
maximum likelihood fit [Tt'I in Mbc and Ai? is per- 
formed on the selected candidates, which are shown 
in Fig. [1] Each signal probability density function 
(PDF) is described by a sum of two Gaussians. For 
the B1 D^TT^ analysis, all three 5° production modes 
{B;b;, B*B^^ and B^^B^J) are fitted simultaneously. For 
the B^ DfK^ mode, only the B*B* component is 
taken into account. The resolutions for Mbc and AE are 
estimated from the MC and scaled by a common factor 
(one for each variable) left free in the DJtt'^ fit. 

Approximating p^* with p^o in the B* B'^j decay, 
the mean values are parameterized, as shown in Table [H 
as functions of the B'^ and B* masses, which are also left 
free in the B'^ D~Tr~^ fit. The continuum (together 
with possible and B^ background) is modeled with 
an ARGUS function [l3| for Mbc and a linear function for 
AE. A non-parametric two-dimensional PDF, obtained 
from MC with the KEYS method [HI , is used to describe 
the shape of the B^ — > D*~tt'^ background. 



TABLE L Parameterization of Afbc and AE mean values. 



Signal 



Mean of (Mbc, AS) 



B*B* 

b:bo 
B°m 



For the B'^ DfTn'^ mode, the three signal yields 
are expressed as a function of three free parameters, 
B [Bg D^n'^), /b'S'i and fg.go, with the relations 
Nm = NboB {B° i)77r+) /m J2k ef'Sfe where M is one 
of the three BY'Bf' -pair production modes and k runs 
over the D~ modes; the third fraction is defined as 
Jbobo = 1 - fB^Bi - fB'BO- The values of J^k^f^k, 
which are the total branching fractions [l^ weighted 
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FIG. 1: (MbcAE) scatter plots for B" DJtv+ (top) and 
B° — > DJK^ (bottom) candidates. The three boxes in the 
top plot are the ±2.5(t signal regions {B^Bg, B^Bg and 
BgB^, from top to bottom) while those in the bottom 
plot are the ±2.5 a B*B* regions for signal (solid) and for 
Bg — > Dj-K^ background (dashed). 




FIG. 2; (a) Afbc distribution of the B^ tt"*" candidates 

with AS in the signal region [-80, -17] MeV. (&) A^: 

distribution of the B^ — > D~tt^ candidates with Mbc in the 
B'B* signal region [5.41,5.43] GeV/c^ The different fitted 
components are shown with dashed curves for the signal, dot- 
ted curves for the D*~tv'^ background, and dash-dotted 
curves for the continuum, (c) and (d) shows the same distribu- 
tions but using the B*B° signal region {AE £ [-57,9] MeV 
and Mbc G [5.38,5.40] GeV/c^). 



TABLE II: Signal efficiencies, yields (A'^) and significances (S). 



T{5S) mode Y^k^^Bk 


N 


S 




— > Ds tt"^ mode 


161 ± 15 






1.58% 


U5t\t 


21. Oct 


BtB° 


1.58% 




2.7cr 


B°BO 


1.56% 


4.0tt.? 


l.lcr 




DfK"^ mode 






B* B* 


1.12% 


6.7t|* 


3.5 cr 



TABLE III: Relative systematic uncertainties (in %) for 
B (B" ^ D-TV+) and B {B° DfK^). 



by the reconstruction efficiencies, are listed in Table El 

Figure [2] shows the Afbc and AE projections in the 
B*B* and in the B*B^ regions of the data, together with 
the fitted function. In the Mbc distribution, the three 
signal components arc present due to overlap of the sig- 
nal boxes; the peak on the right (middle, left) is due to 
B*B* {B*B°, B°^o) production. Table |TT] presents the 
fitted signal yields as well as the significance defined by 
S = a/2 In (£max/'Co) where £max (-Co) is the value at 
the maximum (with the corresponding yield set to zero) 
of the likelihood function convolved with a Gaussian dis- 
tribution that represents the systematic errors. 

Systematic uncertainties on the branching fractions are 
shown in Table IIIII Those on Jb-b* and Jb'bo are 
mainly due to PDF uncertainties. Those due to the beam 
energy, the momentum calibration and the p*^, ~ p^o 
approximation are propagated as systematics on the B* 
mass and B^ mass. The momentum normalization un- 
certainties are much more important in the latter case 
because the measured energy of the -6° candidate is used 
instead of the beam energy. 



Source 




D-7V+ 






Integrated luminosity 


+1.3 


-1.3 


+1.4 


-1.2 


T{5S) 


+4.8 


-4.4 


+5.0 


-4.4 


fs 


+13.3 


-13.3 


+13.6 


-13.4 


fg, g. 






+4.8 


-4.1 


D~ branching fractions 


+6.6 


-6.1 


+6.8 


-5.9 


Efficiencies (MC stat.) 


+1.2 


-1.2 


+1.5 


-1.3 


Efficiencies {R2, cos^hd) 


+4.8 


-4.8 


+4.8 


-4.8 


TT*, if^ identification 


+5.4 


-5.4 


+5.2 


-5.2 


Track reconstruction 


+4.0 


-4.0 


+4.0 


-4.0 


PDF shapes 


+1.0 


-1.0 


+3.3 


-2.7 


Total 


+17.8 


-17.5 


+19.0 


-18.1 



We measure the branching fraction S (S" — > L*^ 7r+) = 
[3-67lo;33(stat.)l^;4^(syst.) ± 0.49(/^)] x lO^^ where the 
largest systematic uncertainty, due to fs, is quoted sep- 
arately, the fraction /^.s. = (90.ll:i^;^ ± 0.2) % and the 
two fitted masses mgo = (5364.4 ± 1.3 ± 0.7) MeV/c^ 
and TOB. = (5416.4+0.4 + 0.5) MeV/c^. These four 
measurements supersede the previous Belle values 'l^. 
We obtain for the first time values for the two fractions 
/b.so = {7.3tli ± 0.1) % and f^ogo = (2.6^11) %, us- 
ing the correlation (—0.77) between fB*B* and fs'BO- 

Our branching fraction is compatible with the CDF 
result [ilH, and is slightly higher (1.3ct) than B{B° 
£)~7r+) [l2|. The value of fg,g, is significantly larger 
than the theoretical expectation of « 70% ^, |6|]. The 
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mass is compatible with the world average value 
while our value for the B* mass is 2.6 a larger than the 
result from CLEO 21]. The mass difference obtained, 
ms* — 771^0 = 52.0 ± 1.5 MeV/cf, is 4. Oct larger than 
the world average of mg*o — toro |l2| . while heavy-quark 
symmetry predicts equal values 

The distribution of the angle between the momen- 
tum and the beam axis in the CM frame is of theoretical 
interest [1] and is presented in Fig. [3] for the signal events 
in the B*B* region, using the sPlot method [22] ■ A fit to 
a 1 -I- acos^fl^o distribution returns x^/n.d.f. = 8.74/8 

and a ~ — 0.591q :[g. It has been checked that the sig- 
nal efficiency does not depend on this angle. We naively 
expect a = —0.27 by summing over all the possible po- 
larization states. 



0.6 -0.2 0.2 0.6 1 
C0S(9*B.) 

FIG. 3: Fitted distribution of the cosine of the angle between 
the Bs momentum and the beam axis in the CM frame for 
the T(55') B*B; signal. 



In summary, a large B'^ —>■ tt"*" signal is observed 
and six physics parameters are measured: the branching 



D: 



= [3.67lO:i(stat.) 



-1-0.43, 
-0.42 



syst.) 

the fractions of the B'^ pair production 

+3.8 . 



fraction B (B° 
±0.49(/s)]xl0- 

modes at the T(55') energy, f^.g, = (90.ll4:o ± 0.2) %, 
fs^So = (7.31^-3 ±0.1)%, /^OB-o = (2.6l|^)%, and 
the masses niB* = (5416.4 ± 0.4 ± 0.5) MeV/c^, niBO = 
(5364.4 ± 1.3 ±°0. 7) MeV/c^. In addition, evidence 
(3.5 fj) for the B^^DJK^ decay is obtained, lead- 
ing to a measurement B {B° -> DfK^) = [2.4tj g(stat.) 
±0.3(syst.) ±0.3(/,)]xlO-4. 
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For the B° -> DfK^ mode, mean values and reso- 
lutions for B° DfK^ and 5° D^tt^ components 
are calibrated using the results of the D~'it^ fit. 

The four yields (signal, continuum, B^ DJtt^ and 
5° £'*^7r+) are allowed to float, but, due to the very 
small contribution of B^ D*~tt'^ , the ratio between 
the yields of i?° -D*"7r+ and B° D~Tr'^ is fixed from 
a fit to data without kaon identification. 

The fit results are shown in Fig. 2] and Table [Til Sys- 
tematic errors are presented in Table Hill We find 6.7^27 
signal events (3.5 cr), corresponding to B (B^ DfK^) 
= [2.4tJ;^(stat.)±0.3(syst.) ± 0.3(/,.)] x lO'^, using 
the previously fitted value of Jb'B'- the ratio 

B (5,0 ^ DJK^) IB [Bl ^ D-'K^) = (6.5t3;5) %, the 
errors are dominated by the low 5° —* DJK^ statistics. 




-0.3 -0.2 -0.1 



(GeV/c^ A E (GeV) 

FIG. 4: Left: Mbc distribution of candidates 
with AS in the B^Bl signal region. Right: Ai5 distri- 
bution of the B° Df candidates with Mbc in the 
BsB* signal region; the left (right) peak is the 5° Df 
(Bj D~it'^) component. The dashed, dotted and dash- 



dotted curves represent the signal, 
grounds, and continuum, respectively. 
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